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We revisit the theory of magnetoresistance for a system of nanoscopic magnetic granules in metal- 
lic matrix. Using a simple model for the spin dependent perturbation potential of the granules, we 
solve Boltzmann equation for the spin dependent components of the non equilibrium electronic dis- 
tribution function. For typical values of the geometric parameters in granular systems, we find a 
peculiar structure of the distribution function of conduction electrons, which is at variance with 
the two-current model of conduction in inhomogeneous systems. Our treatment explicitly includes 
the effects of dipolar correlations yielding a magnetoresistance ratio which contains, in addition to 
the term proportional to the square of uniform magnetization (/i), a weak anisotropic contribu- 
tion depending on the angle between electric and magnetic fields, and arising from the anisotropic 
character of dipolar interactions. 



INTRODUCTION 



Granular magnetic systems displaying giant magne- 
toresistance (GMR) eflFect ||l|, |^, were put forward as 
an alternative to the previously known magnetic mul- 
tilayers H, Q], due to their practical advantages of eas- 
ier fabrication and higher stability. On the other hand, 
the relevant physics in granular systems presents greater 
challenges than in the case of multilayered systems, for 
which classical |^, semiclassical or quantum ^ solu- 
tions are available. Of course, in all these cases the neg- 
ative GMR has the same origin in that the conduction 
electrons with two spin polarizations flow easier at in- 
creasing alignment of locaHzed magnetic moments. But 
a difficulty results from the fact that, while the mag- 
netic state of a multilayered system is described by only 
few classical variables (the orientations of magnetization 
in each magnetic layer), that of a granular system in- 
volves a statistical ensemble of nanoscopic granule mag- 
netic moments which generally cannot be considered in- 
dependent. Also the magnetotransport in 3D granular 
systems cannot be reduced to any of the paradigmatic 
cases in the layered systems, current-in-plane (CIP) or 
current-perpendicular-to-plane (CPP) |^. Nevertheless, 
currently accepted treatments of magnetotransport in 
inhomogeneous materials (both granular and multilay- 
ered) share a common concept: the two-current model, 
in which the spin up and spin down electron subsystems 
carry current independently 1^, 11, |l^. This model 



reasonably reproduces the main features of GMR in gran- 
ular materials. The most common one is the proportion- 
ality of the magnetoresistance ratio Ap/p to the square 
of uniform magnetization, Ap/p cc m? , as first observed 
by Xiao et al |g| and theoretically explained by assum- 
ing the granules to have a single size and to be uncor- 
related |[l2j. The observed proportionality of Ap to the 
inverse of granule diameter d ||l^, |l3| is accounted in these 
models by assuming that the spin dependent scattering 
in mainly at the surface of the granules. Experiments 
also show deviations from behavior |^, |ll|, ^ [T^ 
which have been attributed to granule size dispersion 



||l2| , HJ, [TtI, and/or correlations between the granular 
moments ||11|, p^ , p^ , p^ , pO[ , most pronounced when the 
applied magnetic field is not too strong compared with 
the intergranule (dipolar) interactions. 

However, despite appearances, there are some issues 
open for discussion in the current accounts of magneto- 
transport in granular materials. 

The treatments referred to above usually consider the 
scattering from individual impurities inside the granules 
to be incoherent, despite the fact that electronic mean- 
free path can be much larger than granule diameter or 
even inter granule distance ||l^. In the present paper we 
extend a model put forward by Pogorelov et al |0 in 
which the scattering was considered coherent from the 
whole granule volume and the geometrical factors aris- 
ing from this coherence where properly taken into ac- 
count, albeit within the context of two-current model. 
We now perform a full treatment of the kinetic Boltz- 
mann equation and find that, within this framework, the 
two-current model does not hold. In the presence of dom- 
inantly forward scattering, spin flip processes become a 
crucial determining factor of the electronic distribution 
function. In this model the l/d dependence of Ap is ob- 
tained without assuming specific surface scattering. 

The effect of magnetic correlations on transport is usu- 
ally treated in the assumption that the resistance is pro- 
portional to (/X]^ • /X2) In, H^, |l8| (where /Xj^ and are 
the magnetic moments of granules) , a phenomenological 
result that can be traced back to Gittleman et al ||2l]| . 
The coherent scattering model of Ref. ||l^ put this result 
on a firmer footing. Treating in this model the inverse 
relaxation time as a spin dependent tensor, determined 
by a squared Born scattering amplitude of conduction 
electrons by magnetic granules, one readily gets a cor- 
relation term (/Xj^ • /Xj) due to coherent spin-dependent 
scatterings by the moments /^^ 2 of two different gran- 
ules. Also this theory determines the proper weighting, 
with distance, of the contribution of this correlation to 
transport. Other known approaches to magnetic correla- 
tions have been either phenomenological ||ll|, or numer- 
ical ii,|ii,|ii. 

In the present work we are able to calculate the mag- 
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netic correlation function analytically for temperatures 
well above the characteristic dipolar energy, believed to 
be the dominant interaction in these systems |Q. These 
correlations are included in the full treatment of the 
Boltzmann kinetic equation. Because dipolar interac- 
tions are not isotropic, we predict a weak dependence of 
Ap on the angle 9h between the electrical and magnetic 
fields. This should not be confused with the anisotropic 
magnetoresistance (AMR) seen in systems with a ferro- 
magnetic percolating cluster or arising from lattice dis- 
tortions due to film substrate stresses |^, p^ . 

In Sec. II we discuss the physical parameters for the 
considered system a nd th e limits for validity of the re- 
lated model. In Sec. Ill A the spin-dependent Boltzmann 



spins of rf-electrons (whole bands), the scattering opera- 
tor can be modeled by the form 



equation for this model is formulated and in Sec. IIIB 



the explicit solution is obtained for the distribution func- 
tion, emphasizing the importance of spin fiip scattering 
in a situation where forward scattering dominates. The 
expressions for resistivity p and magnetoresistance ratio 
p with various observable dependencies (including a 
weak dependence of Ap/p on the angle 9^) are presented 
in section IV . Finally, a general discussion and some com- 



parisons to previous theories and experimental results are 
presented in section 0. 

Partial preliminary results have been already reported 
in Refs. giving a general formulation of the 

model and the expression for GMR in a particular geom- 
etry (parallel electric and magnetic fields). Here we give 
a more detailed presentation of the solution of Boltzmann 
equation, generalize it to any angle between electric and 
magnetic fields and present a full discussion of the impli- 
cations of our results. 



II. DEFINITION OF PHYSICAL SYSTEM AND 
MODEL 

We consider a metaUic system of identical magnetic 
spherical granules of diameter d, randomly embedded 
with volume fraction / into non-magnetic metallic ma- 
trix. This is a reasonable approximation to real granular 
alloys, like ¥eCu, QoAu^ etc., where the magnetic compo- 
nent by transition (T) metal is mostly aggregated in stan- 
dard granules (having sizes only slightly dispersed around 
the mean value and shapes of polyhedra faceted along 
low- index crystalline planes, close to spheres). Conduc- 
tance in such systems is mainly realized by s-like elec- 
trons, shared between noble metal (N) matrix and gran- 
ules, and they are scattered by magnetic d-electrons only 
present in the granules (it should be noted that both 
T- and N-atoms occupy the sites in a common crys- 
taUine lattice). At room temperature, d-electrons are 
distributed in split spin sub-bands and provide almost 
saturated, uniform magnetization within each granule. 
However this magnetization can be randomly oriented in 
different granules. Since Fermi s-electron has in general 
different velocities in N- and T-metals and its spin is 
subjected to some contact interaction with the polarized 



(1) 



Here U and / are the parameters of potential and spin- 
dependent scattering, Taa' is the s-electron spin opera- 
tor. The relevant variable for jth granule, the unit vec- 
tor along its magnetic moment, is considered classical 
and invariable at electron scattering events, since the net 
granule moment p,o typically amounts to ~ lO'^pg and 
its coupHng to the environment should be stronger than 
an energy transfer at electron spin fiip. Thus the model, 
Eq. |l|, implies transitions between different spin channels 
due to spin precession in the field of classical magneti- 
zation within a granule, rather than due to less proba- 
ble spin fiips by individual atomic scatterers. The func- 
tion x(r — Rj) in the simplest approximation is 1 when 
|r — Rj I < To, (ro being the radius of a granule) and zero 
otherwise. It is the Fourier transform of this function 
(see below) that effectively accounts for a distinguished 
role of interface scattering, while it was to be specially 
introduced into the incoherent scattering schemes. The 
scattering potential, Eq. |^, is also easily generalized to 
different sized granules. A similar model, which also con- 
sidered coherent scattering by the granules, was formerly 
proposed by Kim et al ||2J. However, these authors as- 
sumed a dipolar (instead of exchange) coupling between 
granule magnetization and charge carriers spin and did 
not include magnetic correlations between granules. 

In what follows, some important relations will be used 
between the characteristic length scales for this prob- 
lem: the Fermi wavelength Ap = 2i:/k-p, the mean gran- 
ule diameter d, the mean intergranule distance D = 

1 /3 

ij^/Qf) d, and the mean free path £ for conduction elec- 
trons. Namely, we consider them to obey the following 
sequence of inequalities 



£ > Z) > d > Af 



(2) 



which is rather reaHstic for experimental systems (see for 
instance [|o)). 

A particular physical consequence of the relation d ^ 
Af is that scattering is dominated by small angles as 
follows from standard diffraction arguments. We shall see 
that in this situation the spin fiip scattering has increased 
importance in determining the structure of the stationary 
electronic distribution function. 



III. BOLTZMANN KINETIC EQUATION 

A. Spin-dependent distribution functions 

We use the description of non-equilibrium electronic 
state of a granular system, related to the scattering po- 
tential, Eq. ^ and subjected to external electric and 
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magnetic fields, in terms of spin-dependent distribution 
functions /ko- obeying the kinetic Boltzmann equation 
(BE) 



E 



1 



Vk X H • Vk/k.+ (3) 



■ X! (/kV - ha) Wko-,k', 
k',o-' 



= 



where Vk is the conduction electron velocity and Wko-.k'tr' 
its transition probability from kV to ku state per unit 
time. In absence of fields, a trivial steady state solu- 
tion holds /kcr = /k — {exp[/9(ek — £f)]}"^, describing 
the spin degenerate Fermi sphere. The electric field dis- 
torts the Fermi surface (FS) by shifting the Fermi sphere 
and the scattering redistributes electrons, a stationary 
distribution resulting from a balance between these two 
mechanisms. The current density by two spin channels 
is given by 



-e^Vk/ko-, 



(4) 



and in absence of spin flip scattering the up and down 
spin FS are independent. In this case, the a spin FS 
distortion is proportional to To- , the corresponding relax- 
ation time, and the conductivity is therefore proportional 
to + Tj^ , a result used in jl^, |l^ and in all calcu- 
lations in the context of two-current model. Evidently, 
scattering between identical spin states is ineffective in 
relaxing the distortion of the FS if the angle of scatter- 
ing is small. The characteristic transport factor 1 — cosO 
in our case can be estimated as 



1 



COS! 



(Ap/rf)' , 



since the diffraction angle 6* « Ap /rf ^ 1 . 

However, if the spin fiip scattering is present, it con- 
tributes very effectively, without the 1 — cos 9 factor, to 
the relaxation of the difference between the up and down 
spin FS distortions. So, when scattering in mostly in the 
forward direction, the spin fiip scattering forces the up 
and down spin FS distortions to be almost identical. One 
then finds that the transport time is given by 



1 

^tr 



- - 

2 



1 

n 



i.e. the rates, and not the times, must be averaged. We 
also include the effect of correlations between granule 
magnetic moments in a way that the scattering kernel 
of BE involves the connected correlation functions 



C^a(q) 




(5) 



where Cf^ = {fif ■ fJ,° 
thermal average and (• • •)r the average over granule po- 
sitions. Correlations induced by the dipolar interaction 



are not isotropic and depend on the angle between q and 
the external magnetic field. As will be shown below, this 
results in a certain dependence of the magnetoresistance 
on the angle between current and magnetic field. 



B. Solution of kinetic equation 

In presence of fields, we define a usual expansion to 
linearize BE 



/kcr = f] 







>k<T ~ /°(£k - 4>kcr), 



that is 0k(T can be interpreted as the FS deformation. 
Then Eq. H is expressed as 



(6) 



where 0k = ^J2a'^^'^- ^ow BE, Eq. ^, for spatially 
uniform steady state and with neglect of orbital effects 
by magnetic field, can be reduced to an integral equation 
for (jjucT 



evk • E = ^ J df^k't^aa' (k, k')<^k', 



(7) 



The angular integration in Eq. 0is over FS and the kernel 

Wcra' (k,k') is 



-(k,k') = E/ ^PFr.."(k,k")<5(r!k-i^k 



i<.')Saa' 



An 



ppT^cT'Ck, k') 



(8) 



where p-p is the Fermi density of states. The transition 
probability density T„„' (k, k') can be written like in Refs. 
Jl7|, p5|, p6| using the Fermi's golden rule 



r..'(k,k') 



2tt 



{Ka\W\k',a') 



(9) 



where V is the sample volume. 

For the common case of point-Hke scatterers, the vaHd- 
ity of Born approximation only requires that the relevant 
energy scales for perturbation, U and /, be small com- 
pared to the Fermi energy ef But in our case, the 
finite size of scatterers also needs some consideration. In 
a naive view, the strength of perturbation could be rep- 
resented by the scattering rates (see below) 



7& = Y^o^ /^F' 



7? 



^VoI^PF, (10) 



where Vq — nd^/B is the granule volume. For a relevant 
choice of parameters: Vb ~ 10 nm"^, U ^ I ^ 0.3 eV, 
denotes the p-p ^ 10 eV~^nm~^, ep 5 eV, one has a large ratio 



h'^j/sp > 1. However, this estimate ignores the fact 
that the scattering is mostly in forward direction. Taking 
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this into account, the Hfetime of a momentum state is 
estimated as : 

^/r - f{ll + iD/ik^df « /(C/2 + I^)kFd/eF. 

Hence the condition of weak scattering h/r <C £f remains 
vaHd even for values of U and / as large as ~ 1 eV. 

Then BE, Eq. can be rewritten as an operator equa- 
tion in the space of functions defined on FS, 



'1^ = 1^ 



(11) 



where |0e) denotes the driving term and \(j>cr) the FS 
distortion for spin a. 

We use the angular momentum basis \im) and define 
the coefficients 



so that 



0k(T = 0fmcrY™(k). 

Note that since the driving term of BE is evk • E cx 
ei'k-BY'j' (k) , only the £ — l^m — component of |(/)e) 
is non zero. As was already mentioned, correlations in- 
duced by dipolar interactions lead to a dependence of the 
scattering kernel on the angle between the momentum 
transfer vector q = k' — k and the external magnetic field, 
and it is convenient to separate out its (small) anisotropic 
part from the main isotropic one: 



r..'(k, k' 

LJaa' (k, k' 



(12) 



U! 



Let us first restrict consideration to the isotropic scat- 
tering kernel, a;^*^/(q), which depends only on the angle 
between k and k'. The resulting operator in Eq. |ll| 
is diagonal in the spherical harmonic basis with z axis 
along the only distinguished direction of electric field E 
(E-basis) . Note that in this case the quantization axis for 
electron spin operator r can be chosen arbitrary. Then, 



presenting the relevant solution as 4'iq„ = 
4q and 77^0* from 



^10 + <^mo we 



obtain its components 



(13) 



Was 010 +^aaVlO 



0, 



all other components of the FS distortion being zero. We 
should emphasize that, in any case, only the £ = 1 com- 
ponents of 0k(T contribute to the current. The rates ap- 
pearing in Eq. 0| are given by 



(0 = i 



,(■0 



^-j2{io\u;i:i,\io), 



(14) 



or, more explicitly. 



^i5:aa'(10p«|10), 



■E 
E 



as 2 / ^ 



UJ 



E 



4tt 

47r 

dfik' 
47r 



ri,*^,(|k-k'|) (i-cos(0kk')), 

r« (|k-k'|) (c7-CT'cOs(0kk')), 

r« (|k-k'|) C7(l-COs(0kk')), 

r« (|k-k'|) (l-aa'cos(0kk'))- 



Solving Eqs. O, we obtain 



•\/ 47r/3euk£'Ttr'' where 



(i) _i _ (j) UJsaUJas 
'^tr ^ss (7) 



(15) 



is the transport relaxation time in isotropic approxima- 
tion. Introducing i^k into Eq. |^ in the same approxima- 
tion: 



i«Y?(k) = evk-Er«, 



we arrive in a standard way at the Drude resistivity 

Top, 



2 (0' 



As stated in Sec. II and expHcitly shown below, the scat- 
tering probability ro-o-'(k, k') is dominated by small an- 
gles 6'kk'- As a result, the integrals in Eqs. [ijthat involve 
the factor 1 -I- cos 9 are larger than those with 1 — cos 9 by 
a factor of [k-pdf . They appear only in uoaa (and cancel 
in ujsa) so that 



J») -1 



(16) 



Considering the definitions of Eq. Q we can easily con- 
clude that 



(17) 



where the rates rf*2, ^ are defined as 

(7(7 



Si) -1 



PF J (|k-k'|)(l-COs(ekk')). 



Eq. |17| expresses the fact, already discussed in Sec. [II A 



that up and down spin FS have almost the same defor- 
mation and so the relaxation rate of the mean of the two 
deformations is just the mean of corresponding rates. For 
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a typical value kyd ~ 40 the neglected term in Eq. [ij 
turns out to be about 0.5% of the term retained. 

The weak anisotropic term of the scattering kernel, 
w^'^^/, due to the correlation between magnetic moments, 
can be easily included into the present treatment in a 
perturbative way. Then Eq. |ll| reads 



and the solution can be written as 



I0e), 



where j^o- ) is a small anisotropic perturbation. To low- 
est non zero order in it, we have 



E 



(18) 



Projecting out the £ =1 and m = components and us- 
ing the facts that lu)^'^, is diagonal and the only non zero 

component in the unperturbed solution 10^^) is that with 
£ = 1, m = 0, we can solve this equations for the coeffi- 
cients 4>^iQ^ = (1O|0CT°'') = (j^'^iQ + (^"Hw ■ A simplification 
similar to that leading to Eq. |l^ also appfies here and 
we obtain 



(a) 
UJss 



where uji'^ is defined in a similar way to wis in Eq. 
Hence the transport time Ttr is changed by a factor 1 — 
uj^ts /oj^ss compared to r/*^ Eq. ^ and to the order of 
accuracy that we are working in Eq. nil we may write 



-1 -„,W 



(19) 



One might question at this point, whether it is legiti- 
mate to include this last correction while neglecting the 
second term in Eq. |l^. It should be stressed however 
that the two terms in Eq. |l^ are of the same order with 
respect to the small parameter l/(fcF(i). So our theory is 
consistently a lowest non zero order theory in this small 
parameter. On the more practical side, we will see that 
for typical parameter values this correction, arising from 
spin correlations, can in fact be more important than the 
terms neglected in Eq. 11^. 



IV. CALCULATION OF TRANSPORT RATES 

A. Isotropic Kernel 

The principles for the calculation of the scattering rates 
have already been spelled out in Ref. [0. The explicit 
squaring matrix element in Eq. ^ is 



(k,a|M^|k',a') 



qd 



(7(7' 



(20) 



where the function ipi.^) = 3(sin(a;) — a;cos(a;))/a;'^ is the 
structure factor of a sphere, the Fourier image of the 
function x(r-), and W^^, — U5aa' + iTaa' ■ Mj- It is due 
to the presence of factors ip'^{qd/2) in the integrals of Eq. 
[l^ that only values of g < d~^ contribute significantly 
to the scattering. Evidently, Eq. ^ should be averaged 
with respect to the random positions Ri of the granules 
and over the thermal distribution of their moments /x^. 
The latter can be done in the basis with z axis along 
the external magnetic field H (H-basis), where only the 
z component of a magnetic moment has non zero average 

We separate out contribution involving a single gran- 
ule, and write averages of products of magnetic moments 
of diferent granules as products of averages plus con- 
nected correlation functions, Eq. |[ It is then straightfor- 
ward to obtain the explicit formulas for diagonal and non- 
diagonal components of isotropic and anisotropic parts of 
the scattering rates, Eq. fa 



PFr«(g) 



2(77,77/ (/z^) {l + g{q)) 



(21) 



and 



pprW(q) = /^2(^)7|q|(q), 



(22) 



/V^^(y)7|C^(q)- 



Here g{q) is the pair correlation function, already calcu- 
lated in Refs. JTtI , 25, using the excluded volume 
approximation 



-(E< 



-iq-(R,-R,; 



){R} 



^0 Jr>d 



-Smqd) 



(a delta function at g = is neglected because it gives a 
zero contribution to the integrals over q) . The correlation 
functions C||(q) = Cz(q) and C±{q) = Cxin) + Cy{q) 
and the relevant average (/i^) are calculated in the same 
approximation in Appendix B. 

The angular integrations appearing in Eq. ^ can now 
be transformed to integrals over the momentum trans- 
fer q. They can be classified in terms of inverse pow- 
ers of the large parameter kpd (see Appendix A). The 
terms, involving the factor 1 -|- cos 9, are of the order 
{luT^Df li^Fd)"^ whereas the factors involving the fac- 
tor 1 — COS0 oc q^ are smaller by a factor l/(fcFd)^ 
of the order {-fl nj) f I {k^df. 

In the absence of correlations the result for the trans- 
port rate is quite simple (see Eq. O) 
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4/ 



(23) 



where 7^ = lu + lj and the constants a — (9/2) In(fcFd), 
and f3 « 4.172 are defined in Appendix A. Other (non- 
magnetic) scattering mechanisms will give rise to an ad- 
ditive background contribution to the resistivity pb so 
that 



Pb + 



4/ 



This result is identical to the first two terms of Eq. (18) 
in Ref. jljl and it gives a magnetoresistance proportional 
to the square of magnetization 



me 4/3/272 



(24) 



Such proportionality of the magnetoresistance to (Aiz)^, 
at this level of approximation is also obtained in all cal- 
culations in the two-current model |lo|, |l^. Differences 
between these approaches will be discussed in more detail 
in Sec. 0. 



B. Contribution of Correlations 

A calculation of the effect of correlations on transport 
requires expHcit expressions for the correlation functions. 
In this article we consider only the high temperature limit 
of dipolar interactions, k^T > p-q/D'^ when the high tem- 
perature expansion is meaningful. The calculation gives 
(see Appendix B) 



- #rl^2(p|;)V^(gd)P.(cos0, 



'q,hj 



(25) 



where £2(5) = [Cis)/s]^ - [C'{s)]^ and C{s) = coth(s) - 
1/s is the Langevin function. The £2 factor is related to 
the common field effect on the magnetization. A more 
subtle field effect on the transport follows from the fac- 
tor ^2(2;) = (3a;2 — l)/2, depending on the angle 9q^h 
between the scattering vector and the external magnetic 
field. This dependence results from the calculation of the 
correlator C(q) using the H-basis (that where z-axis is 
along H) . It eventually introduces a dependence of the 
transport time on the angle 0h between the current and 
the magnetic field (between E- and H-bases). Returning 
again to the E-basis (with H lying in the xz plane) for 
integration in Eq. || and using the results of Eq. |2^, one 
readily arrives at the following expression for the relevant 



rate ljI 



to 



(a) ^ iPlfpl r .A^o/l 



knT' 



3 VokBT 

X [|n°(0k)p-n"(ek')n°(^k)] 



(26) 



drJk' V' (y )V'('Zrf)^2 (cos 0q,h) X 



The arguments o^^p functions contain q = 2kp sm{0^^' /2) 
were ^kk' is the angle between k and k'. The integrals 
over Qk and Qk' in Eq. |2^ are calculated using the ad- 
dition theorem for spherical harmonics to give 

So, finally we obtain for the transport rate 

-^e 4/ 



Pb 



UeG^ {k-pd)'^ 



a7 



(27) 



;C2{-r^)P2{cos9h) 



SOVokBT "^ksT' 



where again a background contribution, pb, from 
non-magnetic scattering mechanisms (other impurities, 
phonons) is included. The magnetoresistance is then 
modified from Eq. 24 to 



nge-^ (fcpa)* \ 30 J ks-l 

(28) 

where To, the characteristic temperature of dipolar inter- 
actions, is 



To 



j4_ 

Voki 



(29) 



The last term in Eq. |2^ describes the deviation from 
^ {pz)"^ behavior of GMR due to dipolar magnetic cor- 
relations. 



V. DISCUSSION AND CONCLUSIONS 

When the coherence between individual scattering 
events inside a granule is disregarded |jl|, fo|, |l^ the re- 
laxation rate of spin a electrons Ao- = — Aq + crAi 
is proportional to the volume fraction of granules, /, and 
Ai is proportional to the mean normalized magnetic mo- 
ment {pz)- Furthermore, if magnetic scattering at the 
granule surface dominates, Aiis proportional to the sur- 
face to volume ratio of the granules i.e. inversely pro- 
portional to d, the diameter. In the two-current model 
the spin up and down electron subsystems carry current 
independently and the conductivity is given by 



nee' 



(n + n) 



or 



p{h) 



A2 
^0 



A? 



Tlee^ Ao 

As a result the magnetoresistance is 



-Ap: 



A? 



UeC^ Ao 
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From this a proportionahty follows of Ap to (/x^) , \/d 
and /, provided the resistance in zero field is also domi- 
nated by surface scattering at the granules (Aq propor- 
tional to / and 1/d ). 

In our model the entire granule scatters coherently. 
We believe this is reasonable, in view of the fact that the 
electronic mean free path can be larger than the granule 
diameter or even the inter-granular distance This 
corresponds to the homogeneous limit of Camblong et al 
1^. Given the large size of the granules, compared to 
the Fermi wavelength, scattering is predominant in the 
forward direction, concentrated in a cone of angular size 
of order (Ap/d)^ ~ l/(A:F(i)^. As we emphasized, spin 
flip scattering then enforces that the up and down spin 
FS distortions stay in step, and the resistivity is given by 



m 



Hee" 



(At 



The interference between scatterings in different gran- 
ules gives a contribution to Acr proportional to (/i.]^ • /X2) 
which, in the absence of correlations yields a {^z) term 
proportional to Since in our calculation surface scat- 
tering is not specifically distinguished, one might expect 
the scattering cross section to be proportional to the vol- 
ume of the granule ^ . However, as we mentioned 
above, it is decreased by a factor of l/(fcFd)^ by a stan- 
dard diffraction argument. The contribution to the re- 
sistivity carries an extra factor l/{k-pdY arising from the 
transport factor 1 — cosO. As a result the resistivity is 
proportional to \/d. 

Given that our model corresponds to the homogeneous 
limit one could question whether it is capable of predict- 
ing the actual values of magnetoresistance ratio observed 
in experiments. 

By defining a resistivity scale 



{k^dy 



we may rewrite Eq. 03 in the form 



77rTn 



Ap = Ap,„ ^{^izy + ^£2(^^)^2(008^,0 

(30) 

The magnetic field dependent factor (within brackets) is 
zero for zero field and unity for saturating field, where 
correlations no longer contribute, so Ap™ is the max- 
imum value of |ApL The zero field resistance can be 
written as (see eq. 



(31) 



This gives an upper bound for the magnetoresistance ra- 
tio 




P(0) 



a 



(32) 



0.4 0.6 
H(T) 



FIG. 1: Field dependent factor in the magnetoresistance as 
a function of magnetic field at various temperatures. Dashed 
lines - without dipolar correlation correction; full lines - with 
inclusion of dipolar correlation correction at parallel electric 
and magnetic fields (To = 107 K, fio = 1492 and / — 0.2 
0)- 



The ratio l3/a w 0.94/lnfcFd varies between 0.31 0.23 
for kpd in the range of 20 ~ 60. This would imply 
|Ap|//3 < 0.3/, that is considerably smaller than what is 
usually observed below room temperature. The relatively 
small value of 13/ a is, however, a consequence of assum- 
ing a sharp granule interface in the function x(r — Rj), 
which leads to a logarithmic factor in a. To check this, 
we have also calculated with smoother density profiles 
and, for values of the interface thickness as small as 10% 
of the diameter a rather reaHstic /3/q; ~ 0.8 is ob- 
tained. For comparison, Wang and Xiao found maximum 
magnetoresistance ratios of the order of 5% at 300 K and 
15 % at 77 K in a sample of Fe2o Aggo with d — 29 
A ||l^. These authors also found p(0) to be inversely 
proportional to d, which means that the second term in 
Eq. |3| should dominate over pb. 

The effect of magnetic correlations on transport is con- 
tained in the second term of Eq. Our treatment ap- 
plies to dipolar interactions in the high temperature limit 
Po/D^ <C fceT, where D is a typical minimum distance 
between granules. Since D'^ — Vb//, this condition is 
equivalent to T ^ /Tq- Using, as an example, the values 
of the sample mentioned above we estimate Tq « 107 K 
and /To ~ 21 K. This shows that the high temperature 
approximation used to calculate the correlator, Eq. |2^, 
is actually valid down to reasonably low temperatures. 

To illustrate the effect of the correlation term, we com- 
pare in Fig. |l| the factor (^0^ + (77rro/30T)/:2(Mo V^bT) 
for Oh — (strongest correlation effect) with {pz^ alone, 
as functions of magnetic field for different temperatures 
(we used the parameters of the sample mentioned above, 
Tq = 107 K and ^0 = 1492 ^b)- As expected, the correla- 
tions give a larger correction (as much as 10%) at lower 
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FIG. 2: Field dependent factor of magnetoresistance ratio, 
as a function of {fJ.z)'^, for T = 40 K and the same sample 
parameters as in Fig. |l[ The effect of dipolar correlations 
is seen in the curvature, which has opposite signs for parallel 
{9h = ) and perpendicular {6h — it/2) electric and magnetic 
fields. Notice that for the angle Oh = arccos l/VS « 55° the 
dipolar correlation correction turns zero and the straight line 
is recovered. 



field. However, at temperatures lower than /Tq, dipolar 
correlations cannot be treated within high temperature 
expansion, and it may happen then, that the zero field 
correlation effect remains finite. This in fact could rise 
the maximum magnetoresistance above the limit of Eq. 

Apm- The anisotropy we predict may already have 
been noticed in CuAg films by Stearns and Cheng [Q. 
It should not be confused with usual AMR, due to skew 
magnetic scattering, observed in more concentrated sys- 
tems [|||. 

To conclude, we would like to stress once more that 
these results follow from a peculiar structure of the elec- 
tron distribution function, related to the fact that mag- 
netic scattering is mostly in the forward direction. The 
relaxation rate of the difference between the up and down 
spin FS is the fastest one, leading to identical up and 
down spin FS distortions. The resulting model differs in 
some ways from the two-current model of magnetotrans- 
port, but seems to account equally well for the basic fea- 
tures of experiments. Some of the ideas presented here 
may also apply in other systems whenever large struc- 
tures, but still smaller than electronic mean free path, 
dominate magnetic scattering. 
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temperatures. A plot of magnetoresistance ratio Ap/ p vs 
{pz)"^ should be a straight line if the correlations are ne- 
glected. In Fig. plwe illustrate this effect by representing 
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due to the anisotropic nature of the dipolar interactions. 
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for the moment components along the direction joining 
the granules but antiferromagnetic for the perpendicu- 
lar components. The correlation correction vanishes and 
linearity of Ap/ p vs {pz)'^ is restored for the particular 
angle 9h = arccos l/-\/3 « 55°. These features, together 
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experimental detection of the dipolar correlation effect 
in GMR of granular metallic systems. It should be men- 
tioned that granule size dispersion also leads to a non- 
linear Ap vs {pz)'^ 



|17| , |19| but always convex and 
6';i-independent. Both these kinds of non-linearity can 
coexist in real granular systems. At last we note that 
the correlation term by Eq. SO is zero at zero magnetic 
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VII. APPENDIX A 

The angular integrations, over the direction of k', 
in Eqs. |l4| can be transformed into integrations over 
the momentum transfer q = 2fcF sin6'kk'/2. The term 
cos6'kk' = 1 — 2(g/2fcF)^ and the solid angle integration 
element, dfik' = 8nqdq/{2kF)^ . Using of Eqs. ^in Eqs. 
nj involves the following integrals: 



HQ) 



Q 



x-^Tp {x)dx = 

= ^[ln2Q + 7E-l-Ci(2Q) + 
sin 2Q 1 — cos 2Q 



Q 



2Q2 



J{Q) 



x^'ip'^{x)ip{2x)dx 



27 
20 



2 In 2 - 1 + 2Ci (2Q) - 2Ci (4Q) 



sin2Q(l-cos2g) 

Q 



3 

4Q2 



1 
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5 + 4cos2(3 — cos4Q cos2Q — cos4Q 
8Q2 + 



where 7e ~ 0.5772 is Euler's gamma. The transition 
rates, Eq. |l^, are simply expressed through these inte- 
grals: 



4/ 



'tr 



7^/ (fcpd) - J{kFd)8f-ff X 



and in the limit of kpd 3> 1 (taking into account the 
asymptotics of integral cosine Ci(x) sinx/x) we have 
the result of Eq. |^ with a — {9/2)\n{kpd) and /? = 
(54/5)(21n2-l)w 4.172. 



VIII. APPENDIX B 

The Hamiltonian for magnetic interactions of the gran- 
ules is 



and the moment-moment correlation function Q 



C,;. 



dC{s] 



ds 

1-3 cos^ 6 



with s = fioh/kBT. These expressions must still be av- 
eraged over the positions of the granules. We assume a 
uniform distribution with the excluded volume constraint 
||17| and obtain for (^x): 



which can be written in terms of the classical demagne- 
tizing factor Nz of the granular sample as 



where h is the external magnetic field appHed to the sys- 
tem (with z axis rotated compared to the geometry in 



Sec. [V B| ) and 11°'^'^ = {Sap — 3ufju'^j)/rfj (with the unit 



vector Uij — Yij/rij) is the dipolar interaction tensor. 
We made a standard high temperature expansion, with 
respect to the dipolar interaction energy, vaHd for tem- 
peratures for which the dipolar interaction is a small per- 
turbation to the non-interacting Hamiltonean. 

To the first order, we obtain for the magnetic moment 



knT 



Finally, the Fourier transform (restricted by the excluded 
volume) C(q) = C||(q) + C±{q) of the correlation func- 
tion is 



C(q) = 



d^rC(r)e 



2q-r 



^0 Jr>d 

y/^£2(s)^(2<zro)P2(cos0q,h) 
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